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Abstract. We address here the problem of extending the Pesin relation among 
positive Lyapunov exponents and the Kolmogorov-Sinai entropy to the case of 
dynamical systems exhibiting subexponential instabilities. By using a recent 
rigorous result due to Zweimuller, we show that the usual Pesin relation can be 
extended straightforwardly for weakly chaotic one-dimensional systems of the Pomeau- 
Manneville type, provided one introduces a convenient subexponential generalization 
of the Kolmogorov-Sinai entropy. We show, furthermore, that Zweimuller's result 
provides an efficient prescription for the evaluation of the algorithm complexity for 
such systems. Our results are confirmed by exhaustive numerical simulations. We also 
point out and correct a misleading extension of the Pesin relation based on the Krengel 
entropy that has appeared recently in the literature. 
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1. Introduction 

One-dimensional chaotic dynamics are usualy characterized by the existence of a 
positive Lyapunov exponent, which indicates exponential separation of initially nearby 
trajectories [Tj. In recent years, we have witnessed a rapid development in the study 
and characterization of dynamical unpredictable systems in which the separation of 
trajectories is weaker than exponential [2]. For these systems, generically dubbed weakly 
chaotic in the physical literature, the conventional Lyapunov exponent vanishes and new 
concepts and ideas for the characterization of dynamical instabilities are necessary for a 
deeper understanding of their global dynamics. Many results of infinite ergodic theory 
[3] come out as powerful tools in this context, shedding light on several apparently 
unrelated results in the physical literature. Among them, the Aaronson-Darling-Kac 
(ADK) theorem [3] has certainly a central role in these problems. 

The main purpose of this work is to extend the well-known Pesin relation j3] 
for the case of weakly chaotic one- dimensional systems, a matter that has attracted 
considerable attention recently (see [5] for references) and even some controversy [6]. 
For usual one-dimensional chaotic systems, the Pesin relation is given simply by h — A, 
with h and A standing, respectively, for the Kolmogorov-Sinai (KS) entropy and the 
usual Lyapunov exponent. We will show that adequate subexponential generalizations 
of the KS entropy and of the Lyapunov exponent will obey exactly the same Pesin-type 
relation, for almost all trajectories. It is important to stress that the existence of such 
generalization is far from intuitive since we are dealing with nonergodic systems for 
which the typical dynamical quantities do depend on the specific trajectory. We also 
show that the extension based on Krengel entropy proposed in [5] for weakly chaotic 
systems is incorrect. The source of the problem is identified and the correct expression 
is presented. We close by comparing our proposed Pesin-type relation based on the 
subexponential KS entropy and the proposal involving the Krengel entropy. 

2. Pesin-type relation and statistics 



We will consider here the general class of Pomeau-Manneville (PM) [TJ maps xt+i = 
f(x t ), with / : [0, 1] ->• [0, 1] such that 



for x — > 0, with a > and z > 1. From the physical point of view, the original 
PM system (z = 2) is paradigmatic since it corresponds to certain Poincare sections 
related to the Lorenz attractor [TJ. Systems of the type ([1]) exhibit exactly the 
kind of subexponential instability for nearby trajectories that we are concerned here: 
5x t ~ Sxo exp(A Q t a ), with < a < 1. A distinctive characteristic of such class of maps is 
the presence of an indifferent (neutral) fixed point at x = 0, i.e., /(0) = and f'{0) = 1. 
The global form of / is irrelevant for our purposes, provided it respects the axioms of 
an AFN-system [Sj. These systems are known to have invariant densities u(x) such that 




(1) 




a = (z — 1) 



i 



(2) 
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near the indifferent fixed point x = [9j. Clearly, the corresponding invariant measures 
diverge for z > 2. In these cases, we have pictorially two qualitative distinct dynamical 
behavior coexisting, namely a laminar regime near x = and a turbulent one elsewhere, 
resulting eventually in nonergodicty and subexponential separation of initially close 
trajectories. It is noteworthy here that it was recently shown that subexponencial 
instability does imply infinite invariant measure [10]. On the other hand, 1 < z < 2 
leads to a finite invariant measure, which is naturally related to ergodicity and positivity 
of the ordinary Lyapunov exponent. 

For intermittent systems like (OQ), the statistics of a given observable $ for randomly 
distributed initial conditions has some peculiar properties. For ergodic systems, the time 
average t -1 J2k=h ^(f k { x )) converges to the spatial average / fid/j,, with dfi = u(x)dx. 
On the other hand, if the system has a diverging invariant measure, the time average will 
typically depend on the chosen trajectory. Nevertheless, the ADK theorem [3] ensures 
in this case that a suitable time-weighted average does converge in distribution terms 
towards a Mittag-Lefner distribution with unit first moment. More specifically, for a 
positive integrable function and a random variable x with an absolutely continuous 
measure with respect to the Lebesgue measure on the interval of interest, there is a 
sequence {a t } for which 



for t oo, where £ Q is a non- negative Mittag-Lefner random variable with index 
a E (0, 1] and unit expected value. Notice that for 1 < z < 2 (the ergodic regime), 
at ~ t, and the corresponding a = 1 Mittag-Leffler distribution reduces to a Dirac 
5-function. For the subexponential regime (z > 2), we have at ~ t a [8] and, by choosing 
d = In | f'\, the ADK theorem assures the convergence in distribution terms towards a 
Mittag-Leffler distribution of the subexponential finite-time Lyapunov exponent 



for t — > oo. The generalized Lyapunov exponent (j3J) plays for intermittent systems the 
same role did by the usual exponent (corresponding to a = 1 in (j3|)) for one-dimensional 
chaotic systems, see [TTJ and references therein for a recent discussion. 

In order to investigate the connection between subexponential instability and the 
corresponding degree of randomness of an intermittent dynamical system like (JTJ, we 
will consider the Kolmogorov-Chaitin concept of complexity pp. Let us assume that 
the phase space of the map (JTJ is partitioned and completely covered by a set of non 
overlapping ordered cells. A given trajectory {x t } generated by the map ([1]) can be 
represented by a sequence of symbols {s t }, which we assume to be integers such that 
St corresponds to the cell where x t belongs. The next step in the analysis consists 
in eliminating redundancies that may appear in {st} by performing a compression of 
information. This can be done, for instance, by introducing the so-called algorithmic 
complexity function C t ({s t }), which is defined as the length of the shortest possible 




(3) 



1 



t-i 
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program able to reconstruct the sequence {s t } on a universal Turing machine PQ. 
Systems that exhibit some degree of regularity are able to generate sequences of symbols 
at a rate higher than needed for recording their programs. For example, a periodic 
sequence can be recreated by replaying the periodic pattern over the total lenght. 
Typically, for these cases, one has Ct ~ Int. On the other hand, if the trajectory 
is completely random, there is no way of reproducing it other than memorizing the 
whole trajectory, resulting in a sequence length that increases linearly in time, i.e., 
Ct ~ t. The finite time KS entropy is defined simply as ht = Ct/t. An important recent 
rigorous result due to Zweimuller [12] unveils the relation between KS entropy and the 
Lyapunov exponent for systems exhibiting subexponential instability. According to this 
result, we have, for almost all initial conditions, 

Ct - W, (5) 



E£i *(/*(*)) ' JW 

for t — > oo, for any observable function where h^f) stands for the Krengel entropy 
[T3] . which can be expressed by the so-called Rohlin's formula [2] 

/»„(/) = / \n\f\dv (6) 
By choosing again $ = In |/'|, we get from ([5]) the surprisingly simple relation 

h { t a) ^\ { t a \ (7) 
for t — > oo and for almost all initial conditions, where 

h<r ] = % (8) 

is the subexponential generalization of the finite-time KS entropy. The relation ([7]) is 
the most natural generalization of the Pesin relation for systems of the type ([1]). From 
the ADK theorem and ©, we have that both h[°^ and converge in distribution 
terms toward the same Mittag-Leflier distribution. Nevertheless, Zweimuller's result is 
indeed stronger, assuring that, for almost all trajectories, h[ coincides with in the 
limit t — > oo. In this way, the relation §5§ does provide an efficient prescription for 
evaluating the algorithmic complexity of a given trajectory for one-dimensional maps, 
namely 

a^£ln|f (/ fc (*))|, (9) 

k=0 

for large t. The power of the prescription ([9]) resides in the fact that it does provide, 
for the systems in question, a computable way for the calculation of the algorithmic 
complexity function Ct, a well-known non-computable function in general [lj. It is 
important also to remind that, for dynamical systems with infinite invariant measure, the 
invariant density, and consequently, the invariant measure, is defined up to an arbitrary 
multiplicative positive constant. In other words, the transformation u> — > £uj (implying, 
in this way, that b — > £6 in ([2])), with £ > 0, does not have any dynamical implication. 
Zweimuller's construction, based in ([5]), is clear invariant under such transformation. Of 
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course, such "symmetry" is broken in the usual ergodic case due to the normalization 
of the invariant measure. 

We notice that the relation ([71) is compatible with the pioneering work of Gaspard 
and Wang [2], where the standard PM map /(x) = x+ax z (mod 1) was considered. They 
argue, in particular, that the algorithmic complexity Ct for the PM map is proportional 
to N t , the number of entrances into a given phase space cell during t iterations of the 
PM map. By invoking some results from renewal theory [15J, one has 

Prob (N t > c^j -+G a (u), (10) 

for < a < 1 and t — > oo, where c is a positive constant and G a stands for the 
one-sided Levy cumulative distribution function with index a. With the change of 
variable u = r£~ 1//Q , where r a = aT(a), we have that the normalized random variable 
£ = Nt/ (Nt) tends toward a Mittag-Leffler random variable with index a and unit first 
moment for t — > oo (see [16] for more details on the relations between one-sided Levy 
and Mittag-Leffler distributions), in perfect agreement with the predictions of the ADK 
theorem. The possibility of estimating the algorithmic complexity function Ct from N t 
also for generic systems of the type (JT)) is indeed confirmed by our exhaustive numerical 
explorations. (See, also, [8] and references therein.) 

3. Numerical simulations 

The ADK convergence of the generalized Lyapunov exponent (jlj) was exhaustively 
checked and confirmed by numerical simulations for different maps in [11]. The 
Zweimuller prescription for calculating the algorithmic complexity ((9]) assures also an 
ADK-like convergence for the generalized KS entropy dHJ). A possible way of testing 
the consistence of the overall picture is to compare the Zweimuller prescription (J9j) 
with other independent prescription for calculating the algorithmic complexity Ct- For 
this purpose, we consider some simple realizations of the general maps of the type (JTJ), 
namely the standard PM case discussed in [2], the Thaler map [9] 

-V(2-2) 

(11) 

mod 1, and the so-called modified Bernoulli map (see [17] for references ) 

f x + 2 z ~ 1 x z , < x < -, 

m = \ i 2 (12) 

[ x - 2 Z ~ 1 (1 - x) z , - < x < 1. 

The modified Bernoulli map f[T2l has indeed two neutral fixed points at x = and x = 1, 
but this does not alter our analysis since we still have a t ~ t a for z > 2 in this case [8]. 
Motivated by the construction introduced in [2], let us consider the standard partition 
of the interval [0, 1] into two cells, A = [0, x*] and A\ = (x*, 1], where x* is the point 
of discontinuity of the maps, i.e., lim^^- f(x) = 1, with < x* < 1. For the modified 
Bernoulli map ( 1121 . one has simply x* = 1/2, while for the other cases the value of x* 
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Figure 1. The cells Aq and A\ for the Bernoulli map (fT2]) with z = 5/2. Notice 
that the dynamics are regular inside each of the cells, with the trajectories departing 
monotonically from the respective fixed points. Nevertheless, the transition for one 
cell to the other is chaotic. The situation is similar for the PM map with a = 1 and 
for the Thaler map (fTTj) , even though for theses cases the partitions are not symmetric 
as the present case. (See also [5]). 



does depend on the map details, in particular on the value of z. The trajectories inside 
both cells A and A\ are typically regular, the turbulent behavior is associated with 
the transition from one cell to the other, see [2] and Fig [TJ Let N t be the number of 
entrances of a given trajectory into the cell Ai during t iterations of the map. Since the 
contributions for Ct arising from the laminar parts of the trajectories contained inside 
the cells are subdominant for large t, it is natural to expect that, for weakly chaotic 
regimes, C% = jNt for large t, where 7 is some proportionality constant, independent of 
the specific trajectory considered, implying, in particular, that C t /{C t ) = N t /(N t ) for 
large t. We could check by numerical simulations that the subexponential KS entropy 
OH]) calculated directly from Nt, namely 

hl a) N , s 

(13) 



(hl a) ) (N t )' 

does converge toward a Mittag-Leffler distribution with unit expected value. Fig. [2] 
depicts a typical example. This convergence is robust and typically fast, see [UJ for 
a recent discussion. Interestingly, Eq. (fl3"j) and its convergence issues do not depend 
on the specific partition introduced above, even though the specific value of 7 does. 
Fig. |3] depicts the relation between the algorithm complexity calculated by using the 
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Figure 2. Distribution of finite time Kolmogorov-Sinal entropy h t calculated from 
(fHj) for the Bernoulli map (TT2]), with z = 28/13 (a = 13/15), for t = 6 x 10 4 and 
2.5 x 10 5 uniformly distributed initial conditions. The histogram was built directly 
from the numerical data, while the dotted line is the corresponding Mittag-Leffler 
probability density computed with the algorithm of [TB]- For convergence details, see 
[TT] , Similar results hold also for the other considered maps. 

Zweimiiller's prescription fl9]) and the values of N t for different partitions. As one can 
see, both quantities are indeed proportional, with very good accuracy, irrespective of 
the partition employed. We also notice that the value of 7 depends on the details of the 
maps, specifically on the value of z and, consequently, of a, see Fig. HI 

4. Final remarks 

We close with some final remarks on the work [5] and what has led to its misleading 
conclusion that 

/»„(/) = a (A<">) (14) 

for systems of the type (CQ). The first observation is that (114p is incompatible with the 
re-scaling transformation u — > which must be a symmetry for the dynamics in the 
infinite measure case. The dynamical quantity on the right-handed side, whatever way 
the average is taken, must be invariant under such transformation, while the left-handed 
side is certainly not, see (jSJ). Such problem can be elucidated recalling that the ADK 
theorem gives (see, for instance, |llj ) 




Figure 3. Graphics of the algorithm complexity C*, calculated by the Zweimiiller 
prescription ([§]), as a function of N t , the number of entrances of a given trajectory 
into the cell A\ = 1], during t — 10 6 iterations of the Bernoulli map (fT2|) with 
z = 28/13. For sake of clarity, only 2500 points is shown for each case (a)-(d), which 
correspond, respectively, to x* — 1/2, 5/8, 3/4, and 7/8. The linear relation is evident. 
The situation for the other considered maps is similar. 



from which Rohlin's formula (E]) for the Krengel entropy implies immediately that 

b \aj sin(Tra) \ /adk ' v ; 

which is the correct relation between Krengel entropy and a dynamically meaningful 
average of subexponential Lyapunov exponents for maps of the type (jTJ. The ADK 
average is not only a dynamically meaningful average, it is essentially the dynamically 
meaningful average for these systems. For instance, the average of the subexponential 
Lyapunov exponents (j3J) calculated for randomly chosen (with any absolutely continuous 
measure with respect to the usual Lebesgue measure on the interval [0, 1]) initial 
conditions x will converge to the ADK average for large t, see [11] for some recent 
applications of this important fact. Notice also that both sides of (|T6|) are invariant 
under the symmetry lo — > £cu. 

A closer inspection of [5] (see, in particular, their Eq. (10)) shows that they, when 
dealing with the continuous time stochastic linear model proposed in [18] , tacitly choose 
a value for £ such that 
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Figure 4. The proportionality constant 7 between Ct and Nt , calculated with respect 
to the standard partition, as a function of a for the PM (a), Thaler (b), and Bernoulli 
(c) maps. For the three cases, the typical uncertainty in 7 is about 1% for samples of 
10 4 trajectories (computed up to t = 10 6 ). The curves are calculated with increments 
of 1CP 2 in the values of a. 

breaking the measure re-scaling symmetry of ( fl6l) and rendering it in its ^-dependent 
form (fl4"l) . However, one could have chosen any other value for £, leading to a distinct 
value of b and to a completely different "relation" between the Krengel entropy and the 
ADK average. Since these relations do depend on some specific multiplicative constant 
of the infinite invariant measure, they have no dynamical meaning. It is interesting to 
notice that N t is also considered as a Mittag-Leffler random variable in [5] by using 
renewal theory in a different manner, but its relation to C t is not stated. Instead, it 
is used as hypothesis that Z)l=o hi \f'{f k {x))\ oc N t in order to conclude that X[ a ^ is 
Mittag-Leffler distributed. Such assumption presumes the convergence 

(A| a) ) ~* (N t ) (18) 

for almost all trajectories, which is indeed correct, but it is a very strong assumption 
without a prior knowledge of Zweimiiller's relation (jSJ). 

We close by noticing that, comparing (J7J) and f fl6|) . it is clear that the subexponential 
KS entropy is the appropriate entropy for extending Pesin relation for weakly chaotic 
systems. Relation ([7]) is simpler than ( I16p and, mainly, it is more powerful since it holds 
for almost all single trajectories, in contrast with (Tl6|) , where a statistic description 
involving many trajectories is necessary (and, moreover, an invariant measure, which 
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usually is not explicitly known, is required for the calculation of Krengel entropy). 
Furthermore, the ergodic transition a — > 1 in ffl6|) is rather awkward in comparison 
with the same transition for the relation ((7j), which is straightforward and natural since 
the Mittag-Leffler distribution tends to a Dirac 5-function for a — > 1. 

This work was supported by the Brazilian agencies CNPq and FAPESP. 
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